The relevance of polarized bZ production at LHC 
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Abstract 

We consider the Z polarization asymmetry Az = (cf(Zr) — a(Zj)j) / \g(Zr) + cf(Zl)) 
in the process of associated bZ production at the LHC. We show that in the Standard 
Model (SM) this quantity is essentially given by its Born approximation, remaining al- 
most unaffected by QCD scales and parton distribution functions variations as well as by 
electroweak corrections. The theoretical quantity that appears in Az is the same that pro- 
vides the LEP1 Z bb forward-backward asymmetry, the only measured observable still 
in some contradiction with the SM prediction. In this sense, Az would provide the possi- 
bility of an independent verification of the possible SM discrepancy, which could reach, if 
consistency with LEP1 measurements is imposed, values of the relative ten percent size. 



1. Introduction 



The Standard Model is confirmed up to per-mille precision by collider data [lj; moreover, very 
recently, Higgs boson signals (2j[3] seems to rise in a narrow mass window around 125 GeV, 
consistently with predictions based on global fits to electroweak data [T]. 

The question arises of whether all the theoretical SM predictions have been confirmed by 
the related experimental measurements. The answer to this fundamental question is nowadays 
that at least one experimental result still appears in some sizable contradiction (roughly, at 3 a 
level) with the SM , i.e. the measurement of the forward-backward asymmetry of bb production 
at the Z peak [I], A b FB . 

In fact, a number of models have been proposed that might cure the discrepancy (see [5H7j 
and references therein). In particular, a slightly embarrassing fact for Supersymmetric models 
is the difficulty that the simplest MSSM version would face to eliminate the discrepancy, as 
exhaustively discussed in Ref. [8]. 

The aim of this paper is that of showing that a specific observable can be defined at LHC 
that would provide essentially a re-measurement of the same LEP1 A b FB quantity, in spite of 
the total difference of the produced final state. This quantity is defined in the production of a 
bZ pair as the ratio of the difference of production cross sections with different (left, right) Z 
polarizations (Z L) Z R ) divided by the corresponding sum. 

We shall first show in section 2 that this quantity is straightforwardly proportional to A b FB 
at the simplest partonic Born level, providing a possible ten percent deviation from its SM 
prediction if the relevant parameters are chosen to reproduce the experimental LEP1 result for 
the asymmetry. In Section 3 we shall derive the special property of our considered quantity, 

1. e. the fact that it remains unaffected, at realistic levels, by variations of the strong scales and 
of the adopted parton distribution functions (pdfs) as well as by electroweak corrections. This 
would represent, in our opinion, a strong motivation to perform an accurate measurement of 
the asymmetry at LHC in a not far future, as qualitatively discussed in the final conclusions. 

2. The Z polarization asymmetry at tree— level 

We shall consider the process of associated production of a Z boson and a single b-quark, 
represented in Figure (H defined at parton level by subprocesses bg —> bZ involving two Born 
diagrams with bottom quark exchange in the s-channel and in the ^-channel. 

This process has been calculated at next-to-leading order in QCD in a previous paper [9] where 
the theoretical uncertainties assessment on cross section calculation have been addressed as well. 

For our purposes we need, though, a derivation of the expressions of the polarized cross 
sections. This requires a number of formulae that we shall briefly show in what follows, starting 
from the calculation of the various quantities performed at the Born level. 

The interaction vertices involved in the diagrams of Figure [U are defined as follow 
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Therefore, the Born invariant amplitude is given by 
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Figure 1: Born diagrams for associated production of a Z boson and a single b-quark. 
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where e, X 1 are the gluon polarization vector and colour matrix, e is the Z polarization vector 
and q = Pb + P g = Pz + Vii s = <fi <l' = Pb ~ Pg = Pb ~ Pz, u = q 2 with the kinematical 
decompositions 

p b = (E b ;0,0,p) , p b = (E' b ;p'sm6,0,p'cos6) , (3) 
p g = (p; 0, 0, -p) , p z = (E z ; -p' sin 9, 0, -p' cos 9) , (4) 

The decomposition of Dirac spinors and polarization vectors leads to 24 helicity amplitudes 
denoted as F\^ T ^ with A = ±|, /i = ±1, r ± ~, // = ±1, referring to 6, g, 6 7 , Z respectively. 

However, in order to explore quickly by hand the properties of this subprocess we can neglect 
rrib/Mz and m^j y/s terms and consider only the following eight non vanishing amplitudes: six 
transverse ones 

2eg 8 g* h y/p 2eg £ g|^/g 

^++++ = e , r = - Q , (7) 

cos I cos | 

o cos f r COS - , x 

F + - + - = 2eg a g* b -=l , F_ + _ + = 2eg s g L Zb ^ , (8) 
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R M Z tan2 f r, L M Z tan2 f 

F + _ ++ = 2e^f 6 cos-.^.-^= 2 - , F_ + __ = 2e^| 6 cos-.-^.-^ , (9) 
and two longitudinal ones 

F + _ +0 = -2^,f b ^.^ , (10) 

VP V s 

z- r sin | M z 
F_ + _ = 2V2eg s g L zb ^.-^ , (11) 



having defined 

^ (12) 



For our analysis it is instructive to consider the Z density matrix 



^ F ^Ti F \flTj ■ ( 13 ) 
XfiT 

A priori there are nine independent Z density matrix elements. However with the above 
Born terms and neglecting again the subleading terms in they reduce to only five ones 

P ++ = 4«|-^M' + ^(^) 2 ) +9 S^) . (14) 

4e 2 9 ,M-^r[/'' + ^('^ l ) 2 )+ 9 S^ 1 ) • (15) 
p M = 8e 2 S 2 ( S I+ S M)sm 2 5(^|) • (16) 
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p+ . = p » + = _ 4e%f ^ M| , p -. = p .- = 4e%M !iD!| M| , (17) 

COS | yp'Sy/s J COS | yp'Sy/S J 

With these powerful but extremely simple mathematical expressions at hand, we can explore 
some physical observables of the process under consideration that keep informations of the Zb 
vertex structure. Let's stick for the moment at the partonic level. The first obvious quantity 
that one can inspect is the subprocess unpolarized angular distribution: with the colour sum 
Eco^ ^ r (ii) = 4, the unpolarized subprocess angular distribution (averaged on gluon and b 
spins and colours) is given by 

d^6 = 7gL^ ^ |i?W ' 2 • (18) 

XfITfI f 

One sees that it is proportional to + p + p 00 ) and, summing the above density matrix 
expressions, solely depends on (#f J 2 + {dzb) 2 '- 
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E l F wl 2 = ((9z b ) 2 + (9z b ) 2 )C di ff , (19) 



with 



In order to separate the gf 6 and g|f 6 contributions, and to check so their possible anomalous 
behaviors, one needs to be sensitive to different density matrix combinations than the sum just 
found in the unpolarized distribution. This can be achieved only keeping track of the final 
Z polarization. The general procedure of its measurement has been described in [nflE] for 
Tevatron processes. The Z polarization can be analyzed by looking at Z decay distributions, 
for example in lepton pairs. It is shown that each density matrix element is associated to 
a specific 9^ (pi dependence. The polarized quantities, therein called a p and a 7 , respectively 
proportional to — p ) and to (p +0 — p _0 ) are the only ones in which the combination 

(dzb) 2 ~ (dzb) 2 appears, as one can check by using the above expressions (11411171) of the density 
matrix elements. They respectively produce lepton angular dependencies of the types cos B\ 
and sin 29 1 cos <\)\ as compared to the unpolarized part proportional to (1 + cos 2 9\). The specific 
generalization of that analysis to the LHC case is under consideration at the moment. 

From this brief discussion, we are naturally led to define the Z boson polarization asymmetry 
A z in bZ production as 



_ a(Z R ) - a(Z L ) _ {g« b y - (g L Zb ? r 

= a(Z R ) + a(Z L ) (g« b y + (g L Zb y P ° l ' { } 



where C vo \ is given as a convolution involving the bottom quark (b and b) and gluon (cj) 
pdfs: 

Cpol — 7 " n\~ ■ (22) 

1 i (Ml tan2|^ 2X 



/3 / 

As one sees, Eq. [21] is simply proportional to the asymmetry parameter At'. 

b ~ (9 L z b ) 2 + (9 R Zb f ' { 6) 

the same quantity that is measured in the forward-backward bb asymmetry in e + e~ annihi- 
lation at the Z pole j4]: 

A» FB = \A e A b , where A e = & ~ & . (24) 

In order to exhibit the relation between A z and Ah without any approximations, we have 
implemented a numerical calculation of the full helicity amplitude retaining the bottom mass 
effects; in our calculation we require a final state b-quark with p T >25 GeV and rapidity 
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Figure 2: Polarization asymmetry A z in bZ production at LHC with y/s =7 TeV. The green 
band displays the ±1 a bounds [4] for the measured asymmetry parameter Ab while the SM 
prediction [4] is shown in red. 



\y\ <2 to reproduce the typical experimental phase space cuts. The gluon and bottom quark 
in the initial state are folded with CTEQ6 p2] parton distribution functions. The polarization 
asymmetry A z in bZ production at LHC with y/s =7 TeV is shown in Figure [2] as function of 
Ab along with the SM prediction [4J (red band) and the measured LEP1 value [4] (green band). 

As can be argued by inspection of Figure El the A z measurement at the LHC could be suffi- 
ciently sensitive to Ab in order to discriminate between LEP1 measurement and SM prediction 
provided that a ~ 8% precision will be achieved on A z measurement at LHC. To better realize 
if such required precision could be reached in the A z calculation we need now to assess the 
effect of its dominant theoretical uncertainties. 

3. Impact of the scale/PDF choices and radiative corrections 

The previous discussion has been performed at the simplest Born level. The next relevant 
question is that of verifying whether the expression of A z remains essentially identical when 
possible sources of theoretical uncertainties or NLO corrections are considered. 

We have proceeded in the following way. First, we have taken into account possible effects 
of either strong scales or pdf variations; as shown in Ref. [9], these variations generate a sensible 
effect, of the almost ten percent relative size, in the total cross section. Next, we have considered 
the possible contribution of NLO electroweak radiative corrections; their effect on the total and 
angular cross section have been determined in a recent paper [13] and found to be possibly 
relevant. 

The dependence of A z on factorization and renormalization scales, /ip and jiR respectively, 
is evaluated by varying their values simultaneously by a conservative factor four with respect to 
the central value; A z is shown in Figure [3] as function of Ab for ji F = fiR — fc/x with /i = M z 
and fc =1, 3 and 1/3. As can be observed from Figure [3] effect of scales variation on A z is below 




5 



1%. However is worth noting that the total cross section dependence on ji R could be strongly 
reduced by using the "Principle of Maximum Conformality" scale-setting (see for instance [H]). 

The asymmetry dependence on the pdf is examined performing the numerical calculation 
with different pdf sets. In Figure HJ we present A z as function of Ah for three different LO pdf 
sets: CTEQ [T2]; MSTW2008 [E] and NNPDF [16]. As one sees, the dependence on the pdf 
set is below 2% while the total cross section can be affected by large variations of order 7% [17] . 

The NLO EW effects on A z deserve a rather different discussion. In principle, these effects 
would not introduce any appreciable theoretical uncertainty, since the values of the involved 
parameters are all known with great accuracy. The goal of their calculation would simply be 
that of offering a more complete theoretical prediction for A z . In fact, it is well known that 
electroweak corrections can have sizable effects on processes involving W or Z production at 
LHC. We have observed it in associated top and W production jT8l[T9] and recent papers on 
W+jet or Z+jet production had also mentioned it, see [13|,[20]. These effects can reach the 
several percent size and even more than ten percent on the subprocess cross sections. This can 
be immediately understood by looking at the simple Sudakov (squared and linear) logarithmic 
terms which affect the amplitudes at high energy [2U[22]. To estimate the size of this type of 
effect at lower energies one also can use the so-called "augmented Sudakov" terms, in which 
constant terms have been added to the logarithmic ones [23]. Using this approach, one can 
immediately be convinced that the polarization asymmetry A z will be essentially not affected 
by these electroweak corrections. Actually, looking at the transverse Born amplitudes, one first 

remarks that because g zh ~ 5 g zb) the dominant amplitudes are and F The other 

ones will contribute to the total cross section by terms suppressed by a factor 1/25. Then, 
applying the Sudakov rules of Ref . p~8lfT9l[2T1423] , one sees that the leading logs associated to 
the bL and Z states are very similar for these two amplitudes. A raw estimate gives effects of 
several percents in the 1 TeV range which should directly affect the cross section. 

However for A Z) dominated by (|F_+_+| 2 — \F | 2 )/(l^-+-+| 2 + W | 2 ) ra tio, the 

common electroweak corrections to each of these amplitudes in the numerator and in the de- 
nominator will cancel out. So a small non zero effect will only come from the smaller amplitudes 
(which contribute by a factor 25 less) and from the small differences due to the subleading (mass 
suppressed) terms. 

Using the augmented Sudakov expressions written in ref. [23] we have checked that the 
effects on A z reach at most the one percent level. In this spirit, we shall consider in this 
preliminary paper the SM NLO electroweak corrections as probably irrelevant. A more complete 
determination of their numerical effect will be given in a forthcoming paper. 

4. Conclusion 

We have shown that the Z polarization asymmetry A z in bZ production at the LHC is strictly 
connected to the well known forward-backward bb asymmetry at Z pole, A b FB) measured at 
LEP1. Our results indicate that A z is almost free from theoretical uncertainties related to 
QCD scale variations as well as to pdf set variations; this property strongly suggests in our 
opinion a measurement of A z at LHC as a unique candidate to possibly clarify the long standing 
puzzle related to the A b FB measurement at LEP1. 

More in general it can be observed that polarization asymmetry observables would be quite 
relevant theoretical observables at LHC, as shown by a recent paper [24], where a polarization 
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Figure 3: Polarization asymmetry A z as function of Ab for three different choice of factorization 
and renormalization scales, respectively \ip and fip = I^r = k/ao with /x = M z and =1, 
3 and 1/3. 




Figure 4: Polarization asymmetry A z as function of Ab for three different choice pdf set as 
described in the text. 
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asymmetry was studied in the context of polarized top production in association with a charged 
Higgs boson, as a possible way of determining the tan/3 parameter in the MSSM. A rather 
general conclusion of our paper is therefore in our opinion that measurements of polarization 
at LHC would represent a tough but possibly quite rewarding experimental effort. A more 
complete discussion of Z polarization measurements at LHC will be treated in a forthcoming 
paper. 

Acknowledgements. We thank M. Cobal and M. Pinamonti for several discussions con- 
cerning the measurement of polarization at LHC. 
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